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2"^  Japan-US  Symposium  on  Finite  Element  Methods  in 
Large  Scale  Computational  Fluid  Dynamics 


Abstract 

A  summary  of  the  2"^  Japan-US  Symposium  on  Finite  Element  Methods  in 
Large  Scale  Computational  Fluid  Dynamics,  conducted  Mar.  14-16,  1994  at  Chuo 
University  in  Tokyo,  Japan  is  presented.  Abstracts  of  all  symposium 
presentations  are  included.  This  report  is  based  upon  information  collected  via 
symposium  attendance,  review  of  the  symposium  proceedings,  and  a  post 
symposium  conversation  with  Professor  Mutsuto  Kawahara,  Department  of 
Civil  Engineering,  Chuo  University,  who  served  as  the  Japanese  co-chair  of 
the  symposium  organizing  committee. 

1.  Symposium  Background,  Program  and  Registration  Statistics 

The  2*^^  Japan-US  Symposium  on  Finite  Element  Methods  in  Large  Scale 
Computational  Fluid  Dynamics  was  conducted  Mar.  14-16,  1994  at  the  Surugadai 
Memorial  Hall  of  Chuo  University  in  Tokyo,  Japan.  As  the  symposium  title 
suggests,  the  event  was  the  second  in  a  series  of  biennial  joint  Japan-US 
CFD/FEM  symposia.  The  first  was  held  at  the  University  of  Minnesota  in  August, 
1992  and  the  third  is  tentatively  planned  for  1-3  April,  1996,  also  at  the 
University  of  Minnesota.  The  symposium  series  is  one  of  the  products  of 
several  years  of  collaborative  activity  between  a  group  of  Japanese 
researchers  centered  around  the  Civil  Engineering  School  at  Chuo  and  the  US 
Army  High  performance  Computing  Center  located  at  the  University  of 
Minnesota. 

The  three  co-chairmen  of  the  2"^  Japan-US  Symposium  on  Finite  Element 
Methods  in  Large  Scale  Computational  Fluid  Dynamics  organizing  committee 
are  (including  contact  information) 

Professor  Mutsuto  Kawahara 

Department  of  Civil  Engineering 

Faculty  of  Science  and  Engineering 

Chuo  University 

1-13-27  Kasuga 

Bunkyo-ku 

Tokyo,  112  Japan 

Tel:  (03)3817-1811/1807 

Fax:  (03)  3817-1803 

e-mail:  kawa@civiLchuo-u.ac.jp 


Professor  Tayfun  Tezduyar 

Professor  of  Aerospace  Engineering  and  Mechanics 
University  of  Minnesota 

Interim  Director,  Army  High  Performance  Computing 

1100  Washington  Ave.  South,  Suite  101 

Minneapolis,  MN  55415 

Tel:  612-626-8095 

Fax:  612-626-1596 

e-mail:  tezduyar@ahpcrc.umn.edu 
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Professor  Thomas  J.  R.  Hughes 
Division  of  Applied  Mathematics 
Stanford  University 
Stanford,  California  94305-4040 
Tel;  (415)  723-2040 
Fax:  (415)  723-1778 

The  symposium  program  was  divided  into  a  total  of  twelve  sessions,  two 
morning  and  two  afternoon  sessions  on  each  of  the  three  days.  With  the 
exception  of  session  4,  all  were  devoted  entirely  to  presentation  of  papers. 
Session  4  (Monday  afternoon,  14  Mar)  consisted  of  a  panel  discussion,  and  is 
described  further  in  section  2  below.  A  total  of  47  papers  was  presented,  of 
which  17  were  delivered  by  US  based  symposium  participants.  The  remaining 
30  were  presented  by  Japanese  participants.  All  presentations  were  delivered 
in  English  and  all  papers  and  presentation  aids  were  prepared  in  English. 

A  total  of  143  registrants  attended  the  symposium.  Of  the  indicated  total,  23 
registrants  are  US  based  and  the  remaining  120  are  Japanese. 

The  symposium  proceedings  is  a  collection  of  extended  abstracts  of  the  47 
paper  presentations.  A  post  symposium  contact  with  Professor  Kawahara 
indicates  that  full  length  versions  of  all  papers  are  scheduled  to  be  published 
in  special  issues  of  two  separate  journals;  (1)  an  upcoming  special  issue  of  the 
"International  Journal  of  Numerical  Methods  in  Fluids",  John  Wiley  and  Sons, 
and  (2)  the  premier  issue  of  the  new  journal  "International  Journal  of 
Computational  Fluid  Dynamics",  Gordon  and  Breach  Science  Publishers. 
According  to  Professor  Kawahara,  publication  dates  for  the  special  issues  are 
not  currently  available. 

2.  Panel  Session 

Session  4  of  the  symposium  was  devoted  to  a  panel  discussion.  Following 
introduction  of  the  three  invited  panelists  and  a  brief  description  of  the  panel 
format,  each  panelist  spoke  for  about  15-20  minutes,  and  then  entertained 
questions  for  a  few  minutes  further.  The  panel  was  composed  of  Dr.  Jagdish 
Chandra,  Director,  Mathematical  and  Computer  Sciences  Division,  US  Army 
Research  Office,  Professor  N.  Satofuka,  Mechanical  Systems  Engineering, 

Kyoto  Institute  of  Technology  and  Dr.  Charles  J.  Holland,  Director,  Mathematics 
and  Computer  Sciences,  US  Air  Force  Office  of  Scientific  Research. 

Dr.  Chandra  spoke  on  the  topic  "High  Performance  Computing  Issues  and 
Challenges".  He  identified  a  list  of  six  major  application  domains  (i.e.  Advanced 
Material  Development  and  Manufacturing,  Computational  Fluid  Dynamics, 
Biotechnology,  Environmental  Sciences,  Signal  and  Image  Proeessing  and 
Modeling  and  Simulation)  and  then  highlighted  Parallel  and  Distributed 
Computing,  Algorithm/Software  Development,  Graphics  and  Visualization  and 
Human  Resources  as  the  key  technical  issues  in  their  effective  pursuit.  Dr. 
Chandra's  concluding  remarks  summarized  a  vision  to  which  he  referred  as 
"Towards  Teraflop  Computing". 

Professor  Satofuka  chose  "Trends  in  High  Performance  Computing  in 
Computational  Fluid  Dynamics"  as  his  topic.  He  noted  that  his  background  is  in 
finite  differences,  and  he  identified  Grid  Generation,  Algorithmic  Structures, 
Computer  Resources  and  Physical  Models  as  the  key  elements  of  High 


2 


Performance  Computing  in  Computational  Fluid  Dynamics.  He  devoted  an 
appreciable  amount  of  his  discussion  to  structured  and  unstructured  grid 
generation  techniques. 

Dr.  Holland  spoke  on  "Critical  Issues  Facing  High  Performance  Computing". 
Following  is  a  synopsis  of  his  remarks  in  the  form  of  quotations  from  his 
personal  notes.  The  quoted  material  is  enclosed  in  brackets. 

["I  would  like  to  make  three  points  that  will  impact  the  use  of  high 
performance  computing,  and  which  will  influence  the  directions  of  computer 
science  research  at  AFOSR. 

"(1)  We  need  to  think  about  developing  more  focused  application  software 
capable  of  use  at  the  departmental  level  in  industry  that  is  flexible  (open) 
enough  to  incorporate  hierarchical  modeling  and  the  latest  numerical 
advances.  A  lot  of  general  purpose  finite  element  codes  do  not  permit  enough 
flexibility  to  permit  correct  modeling  of  key  phenomena.  Examples  include  the 
treatment  of  fasteners,  laminated  composites.  This  limitation  restricts  the 
ability  to  use  these  packages  on  real  problems.  Hierarchical  models  that  permit 
control  of  the  idealization  error  are  needed  as  well  as  methods  that  allow  a 
priori  analysis  of  the  numerical  error.  We  are  supporting  these  ideas  at 
Washington  University  through  the  support  of  B.  Szabo. 

"(2)  Moving  beyond  simulation,  optimal  multidisciplinary  design  is  another 
potential  application  of  high  performance  computing.  In  the  past,  people  have 
attempted  to  simply  graft  together  their  favorite  simulation  package  with  a 
standard  optimization  package.  This  often  leads  to  non-convergence  and 
errors  in  the  calculations!  I  believe  that  we  need  to  bring  together  finite 
element  modelers  and  optimization  experts  at  the  beginning  to  correctly 
design  optimal  design  packages.  We  are  experiencing  some  success  in  this  area 
with  a  multidisciplinary  research  team  at  VPI  involving  Max  Gunzberger  and 
John  Burns.  They  have  worked  on  an  "artificial"  forebody  design  for  use  in 
wind  tunnel  simulation  at  AEDC. 

"(3)  For  progress  in  computational  fluid  dynamics,  we  will  need  techniques  for 
accurate  long  time  integration  of  unsteady  calculations  in  order  to  eompute 
the  finite  details  of  shock-vortex  calculations  and  other  important  time 
dependent  phenomena.  There  are  a  range  of  numerical  basics  that  need  to  be 
treated  including  the  correct  numerical  treatment  of  time  dependent 
boundary  conditions.  I  would  like  to  mention  some  recent  work  of  David 
Gottlieb  and  colleagues  at  Brown  University.  They  have  developed  a  technique 
that  resolves  the  Gibbs'  phenomenon  which  is  a  problem  with  discontinuities 
in  spectral  methods.  In  flow  calculations,  this  results  in  spurious  oscillations 
around  shocks.  He  has  shown  how  to  use  the  spectral  information  in  a  post 
processing  scheme  to  create  a  new  expansion  in  terms  of  Gegenbauer 
polynomials  that  converges  uniformly  permitting  an  accurate  capturing  of 
the  shock  phenomena.] 

3.  Comments 

The  English  language  competence  of  the  Japanese  participants  in  this 
symposium  is  very  noteworthy.  It  significantly  exceeds  the  norm  for  joint 
Japanese-US  conferences.  The  papers  are  well  written  and  without  exception 
the  English  language  delivery  excellent. 
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The  two  central  research  communities  represented  at  this  symposium,  the  Civil 
Engineering  Department  at  Chuo  University  and  the  Army  High  Performance 
Computing  Research  Center  at  the  University  of  Minnesota,  have  developed  a 
very  effective  working  relationship  over  recent  years.  Many  of  the  Japanese 
researchers  have  visited  AHPCRC  on  sabbaticals,  etc.  and  the  U.  S.  researchers 
(including  the  U.  S.  Army  participants)  are  well  aware  of  who  they  are  and 
what  they  do.  However,  there  seems  to  be  little  or  no  participation  by  the  other 
U.  S.  military  service  laboratory  organizations. 

Although  the  symposium  name  indicates  an  international  CFD  scope,  the 
participants  are  more  representative  of  a  collection  of  FEM  researchers  who 
have  moved  from  an  historical  interest  in  structural  mechanics  into  CFD.  The 
community  and  the  symposium  series  might  benefit  from  a  move  to  invite 
participation  by  representatives  from  other  CFD  methods  to  future  events. 

The  US  computer  science  community  is  very  interested  in  tracking  large  scale 
computing  applications  in  Japan,  and  the  symposium  committee  responded  by 
scheduling  post  symposium  visits  to  several  .Japanese  industrial  computing 
facilities.  However,  with  perhaps  one  notable  exception  (a  visit  to  Hitachi), 
very  few  examples  of  large  scale  computing  activity  were  presented.  The 
impression  of  the  participants  is  that  the  work  exists  but  was  not  on  display. 

4.  Abstracts 

The  symposium  proceedings  include  "extended  abstracts"  of  all  47  presented 
papers.  As  noted  in  section  1  above,  the  symposium  organizers  intend  to  have 
full  length  versions  of  all  47  presentations  published,  approximately  one  half 
in  an  upcoming  special  issue  of  the  International  Journal  of  Numerical 
Methods  in  Fluids  and  the  remainder  in  the  premier  issue  of  International 
Journal  of  Computational  Fluid  Dynamics.  According  to  the  symposium 
organizing  committee,  the  required  arrangements  have  been  made  with  both 
journal  editorial  boards.  However,  discussions  during  the  course  of  the 
symposium  indicate  that  many  of  the  full  length  papers  are  still  in 
preparation.  Thus,  it  might  well  turn  out  that  something  less  than  all  full 
length  papers  are  ever  published. 

Following  is  a  compilation  of  selections  from  the  extended  abstracts.  Typically, 
depending  upon  the  organization  and  content  of  the  paper  being  excerpted,  an 
introductory  section  (in  its  entirety)  accompanies  title  and 

author/organization  information.  However,  in  some  cases,  due  to  length  and  or 
depth  of  content,  the  introductory  section  is  abbreviated.  Also,  in  those  cases 
where  the  excerpted  selection  cites  entries  in  an  accompanying  reference  list, 
the  corresponding  reference  list  is  included. 
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Another  Stabilized  QlQl  Incompressible  Na  vier-Stokes  Solver 
Using  Explicit  Time-Marching* 

P.M.  GRESHO  and  S.T.  CHAN 
Lawrence  Livermore  National  Laboratory,  Livermore,  CA,  USA 


ABSTRACT 

While  on  our  way  to  investigating  the  potential  utility  of  a  semi-implicit  time-stepping  scheme  based 
on  a  second-order  projection  method  in  which  the  Q]Q]  element  (bilinear  in  both  velocity  and  pressure)  is 
stabilized  by  the  simple  expediency  of  replacing  the  (unstable)  matrix  representing  the  Laplacian  associated 
with  the  consistent  discrete  projection  by  a  stable  but  less-consistent  matrix  in  which,  for  example, 
spurious  pressure  modes  are  traded  in  for  imperfect  mass  conservation,  we  have  spent  some  time  studying 
and  testing  a  simpler  but  related  notion:  forward  Euler  time-stepping  in  which  the  (same)  consistent 
discrete  matrix  associated  with  the  pressure  Poisson  equation  (PPE)  is  again  replaced  by  the  same  less- 
consistent  matrix.  While  the  final  goal  is  to  test  the  'approximate'  projection  method  in  which  a  second- 
order  Adams-Bashforth  method  is  applied  to  the  advection  terms  and  trapezoid  mle  to  the  remaining  terms, 
and  the  (approximate)  projection  matrix  problem  is  solved  via  the  world's  best  multi-grid  scheme,  this 
paper  reports  on  our  interim  goal  of  seeing  how  cost-effective  the  same  matrix  replacement  would  be  for  the 
simplest  code  to  write,  debug,  and  understand. 

But— as  is  so  often  the  case  with  incompressible  viscous  flow— even  forward  Euler  is  not  as  simple  as 
we  had  assumed;  e.g.  the  Differential  Algebraic  Equations  (DAE's)  associated  with  the  Q]Q]  forward  Euler 
scheme  are  not  nearly  as  'useful'  as  one  might  expect,  and  the  hopes  of  precluding  all  pressure  modes  may 
not  be  realized.  In  fact,  we  have  even  learned  that  the  simple  and  obvious  one-for-one  matrix  replacement 
starting  with  the  original  Index  1  DAE's  does  not  yield  a  method  that  'works'— a  statement  that  is  easy  to 
support-and  we  do  so  below. 

♦  Work  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by  the  Lawrence  Livermore 
National  Laboratory  under  contract  No.  W-7405-Eng-48. 


Analysis  of  Shear  Layers  Based  on  the  Lattice  Gas  Model 
Atsushi  YU  ASA,  Takaji  INAMURO*  and  Takeshi  ADACHI 
Advanced  Technology  Research  Center, 

Mitshbishi  Heavy  Industries,  Ltd.,  8-1,  Sachiura  1-chome,  Kanazawa-ku,  Yokohama  236,  Japan 
1.  INTRODUCTION 

Conventional  numerical  analyses  of  fluids  are  usually  performed  by  using  the  NavierStokes  equation 
derived  from  the  continuum  theory  in  which  fluid  is  macroscopically  regarded  as  the  continuum.  As  an 
opposite  approach,  the  "lattice  gas  model  (LGM)"  based  on  a  microscopic  physical  model  called  "cellular 
automaton"  is  recently  paid  much  attention  to  as  a  new  computational  methodn-2].  The  LGM  imitates 
motion  of  particles  in  a  very  simplified  microscopic  world  with  completely  discrete  nature.  It  has  begun  to 
be  applied  to  numerical  simulations  for  various  complicated  physical  phenomena[3-5]  such  as  fluid 
motion,  molecular  diffusion,  chemical  reaction  and  so  on. 

In  this  paper,  we  used  the  two-color  LGMt-^‘^1  and  evaluated  the  validity  of  this  model  by  applying  it 
to  the  analysis  of  shear  layers. 

*  He  is  currently  working  with  the  Department  of  Chemical  Engineering,  Kyoto  University. 
REFERENCES 

[1]  S.  Wolfram,  J.  Stat.  Phys.,  45  (1986)  471-526. 

[2]  U.  Frisch,  et  al.,  Complex  Systems,  1  (1987)  649-707. 
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[3]  D.  dHumieres,  and  P.  Lalletnand,  Complex  Systems,  1  (1987)  599-632 

[4]  C.  Burges,  and  S.  Zaleski,  Complex  Systems,  1  (1987)  31-50. 

[5]  D.  d'Humieres,  P.  Lallemand,  and  G.  Seaiby,  Complex  Systems,  1  (1987)  633-647. 


Finite  Element  Analysis  for  Three-Dimensional 
High  Reynolds  Number  Flows 

K.  KAKUDA  and  N.  TOSAKA 

Department  of  Mathematical  Engineering, 

College  of  Industrial  Technology, 

Nihon  University,  Narashino,  Chiba  275,  Japan 


1.  INTRODUCTION 

It  is  well  known  that  the  conventional  Galerkin  finite  element  and  centered  finite  difference 
approximations  lead  to  sparrows  oscillatory  solutions  for  flow  problems  at  high  Reynolds  numbers.  To 
overcome  such  oscillations,  various  upwind  schemes  have  been  successfully  presented  in  both 
frameworks[l]. 

In  our  previous  work,  we  presented  a  finite  clement  scheme  based  on  the  PetrovGalerkin  weak 
formulation  using  exponential  test  functions  for  solving  the  two-dimensional  incompressible  Navier-Stokes 
equations  at  high  Reynolds  numbers[2].  The  method  was  tested  on  a  square  cavity  flow  and  flow  around  an 
obstacle,  and  the  numerical  results  were  compared  with  those  of  the  other  existing  ones.  Consequently,  we 
demonstrated  that  the  method  was  capable  of  solving  the  incompressible  NavierStokes  equations  accurately 
and  in  a  stable  manner  for  high  Reynolds  numbers. 

The  purpose  of  this  paper  is  to  extend  the  Petrov-Galerkin  finite  element  method  using  exponential 
test  functions  to  three-dimensional  incompressible  viscous  flow  problems  at  high  Reynolds  numbers.  In 
order  to  show  the  workability  and  the  validity  of  its  extension  to  three-dimensional  problems,  flows  in  a 
cubic  cavity  up  to  Reynolds  number  of  10^  are  demonstrated. 

REFERENCES 

[1]  R.  Peyret  and  T.D.Taylor,  Computational  Methods  for  Fluid  Flow,  SpringerVerlag,  New  York 
(1983). 

[2]  K.Kakuda  and  N.Tosaka,  Finite  element  approach  for  high  Reynolds  number  flows.  Theoretical  and 
Applied  Mechaiucs,  41  (1992),  pp. 223-232. 

[3]  A.K.Prasad,  C.  Y.Pemg  and  J.R.  Koseff,  Some  observations  on  the  influence  of  longitudinal 
vortices  in  a  lid-driven  cavity  flow,  AIAA  paper,  88-3654-CP,  (1988) 


Application  of  FEM/FDM  Overlapping  Scheme  to 
Heat  and  Mass  Transfer  Problem 

M.  IKEGAWA,  M.  KAIHO  and  C.  KATO^ 

1  Mechanical  Engineering  Research  Laboratory,  Hitachi,  Ltd.,  Tsuchiura,  Ibaraki,  Japan 
1.  INTRODUCTION 

There  are  two  major  numerical  analysis  techniques  that  are  commonly  applied  to  computational  fluid 
dynamics:  the  finite  element  method  (FEM)  and  the  finite  difference  method  (FDM).  Finite  element 
method  is  best  suited  for  accurate  simulations  of  flows  within  arbitrarily  shaped  flow  geometry,  while 
FDM  analysis  is  best  suited  for  minimizing  computational  time  and  storage  requirements.  To  effectively 
solve  flow  in  a  complex  region  of  interest  effectively,  Nakahashi  and  Obayashi  presented  an  FDM/FEM 
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zonal  approach  to  analyze  compressible  flows  in  turbine  cascades  [  1  ] .  In  that  approach,  the  region  near  wail 
boundaries  was  covered  with  a  boundary-fit  grid  and  the  remaining  regions  were  covered  with  a  finite- 
element  mesh.  The  authors  propMJsed  a  different  technique  based  on  an  FEM/FDM  overlapping  scheme  for 
viscous,  incompressible  flow,  which  has  been  applied  to  a  two  dimensional  moving  boundary  problem 
around  high  speed  trains  [2],  In  this  technique,  FEM  meshes  are  used  in  the  region  near  wall  boundaries 
and  FDM  meshes  are  in  the  other  regions.  ITus  technique  permits  accurate  simulation  of  boundary  layers 
and  solution  of  moving  boundary  problems,  with  a  minimum  of  computer  resources.  From  a  mesh 
generating  point  of  view,  this  scheme  is  thought  to  be  better  than  conventional  schemes,  because  it  is 
relatively  easier  to  generate  an  FEM  mesh  near  the  complex-shaped  wail  compared  with  generating  an  FEM 
mesh  for  the  entire  calculation  domain.  In  this  paper,  some  numerical  results  for  heat  and  mass  transfer 
obtained  with  this  calculation  scheme  are  presented. 

REFERENCES 

[1]  Nakahashi,  K.  and  Obayashi,  S.,  FDM-FEM  Zonal  Approach  £  or  Viscous  Flow  Computations  Over 
Multiple  Bodies,  AIAA  paper  No.  87-0604,  1987. 

[2]  Ikegawa,  M.  et  al.,  FEM/FDM  Composite  Scheme  {or  Viscous  Incompressible  Flow  Analysis,  to  be 
published  in  Comput.  Methods  Appl.  Mech.  Engrg. 

[3]  VieceHi,  J.A.,  A  Computing  Method  for  Incompressible  Flows  Bounded  by  Walls,  J.  of 
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1.  INTRODUCTION 

We  present  massively  parallel  finite  element  computations  of  compressible  flows  with  the  Connection 
Machine  CM-5  and  CM-200.  The  flow  problems  we  consider  typically  come  from  aerospace  applications. 
We  use  conservation  variables  formulations  of  the  Euler  and  Navier-Stokes  equations,  and  to  prevent 
potential  numerical  oscillations,  we  use  streamline-upwind/Petrov-Galerkin  (SUPG)  stabilization.  For 
compressible  flows,  this  stabilization  technique  was  first  used  by  Tezduyar  and  Hughes  [  1  ].  In  recent 
years  this  technique  was  enhanced  by  a  shock-capturing  term  [2],  and  this  enhancement  rendered  this 
technique  veiy'  comparable  in  accuracy  to  those  that  followed  [1]  and  formulated  in  entropy  variables  [3], 
The  numeric^  strategy  we  use  here  to  solve  compressible  flow  problems  with  moving  boundaries  and 
interfaces— the  deformablespatial-domain/stabilized-space-time  formulation— was  first  introduced  in  the 
context  of  incompressible  flows  [4],  and  soon  after  that  was  applied  to  compressible  flows  [5]. 

Our  massively  parallel  computations  of  compressible  flows  started  in  early  1 992,  with  our  first  results 
reported  in  Tezduyar  et  al.  [6].  A  later  version  of  reference  [6]  appeared  as  a  journal  article  by  Tezduyar  et 
al.  [7].  Our  parallel  implementations  are  based  on  the  assumption  that  the  mesh  is  unstructured,  and  this 
makes  our  parallel  computations  applicable  to  practical  problems  with  complicated  geometries.  For  time- 
integration,  we  typically  use  implicit  methods,  although  we  also  have  the  capability  to  use-and  have  in 
fact  used  in  some  problems — explicit  methods.  Today  we  are  at  a  point  where  we  have  the  capability  to 
solve  a  large  class  of  practical  problems,  including  those  involving  moving  boundaries  and  interfaces  and 
those  in  3D,  using  implicit  formulations  with  problem  sizes  over  5,000,000  equations  and  computational 
speeds  at  12  GigaFlops.  We  believe  this  new  capability  has  pushed  large-scale  finite  element  computation 
of  compressible  flows  to  a  new  era  in  numerical  simulation. 
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1.  INTRODUCTION 

Temam  and  his  colleagues  [1]  have  been  developing  the  dynamical  systems  theory  for  dissipative 
dynamical  systems.  It  is  shown  that  approximate  inertial  manifolds,  which  represent  interaction  between 
large  and  small  scale  motions,  approximate  the  global  attractors  of  such  dynamical  systems.  In  order  to 
investigate  their  own  theory,  named  nonlinear  Galerkin  methods,  they  have  carried  out  computations  using 
a  hierarchical  finite  element  for  the  2D  Burgers  equation.  We  will  show  a  different  finite  element  to 
implement  the  nonlinear  Galerkin  method  for  the  incompressible  Navier-Stokes  equations. 
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Space-Time  Finite  Element  Computation 
of  Compressible  Flow  Between  Moving  Components 

G.P.  WREnI,  S.E.  RAY2,  S.K.  ALIABADI^,  and  T.E.  TEZDUYAR^ 

^U.  S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  USA 
^University  of  Minnesota,  Minneapolis,  MN,  USA 

Recent  advances  in  finite  element  strategies  to  solve  compressible  flow  problems  have  allowed  the 
solution  of  equations  describing  challenging  applied  problems.  The  problems  investigated  in  this  paper 
involve  nearly  incompressible,  barotropic  flow  (i.e.,  flows  in  which  small  changes  in  density  lead  to  large 
changes  in  pressure).  Such  problems  challenge  the  robustness  of  the  numerical  method,  and,  in  some 
cases,  limited  experimental  data  with  which  to  compare  the  solution  are  available.  In  general,  the  physics 
and  dynamics  of  nearly  incompressible  flows  at  high  pressure  over  complicated  geometries  are  not  yet  fully 
understood.  In  particular,  better  understanding  of  thermal  characteristics  and  the  jxttential  of  the  flow  to 
remain  attached  to  the  boundary  and  to  develop  vortices  is  needed  to  design  injectors  for  such  flows. 

The  finite  element  method  is  ideally  suited  for  the  class  of  problems  described  previously  because  of 
the  inherent  geometric  flexibility  of  the  method.  Since  grid  indexing  is  not  necessary,  unstructured  grids 
which  are  often  needed  to  effectively  solve  for  flow  properties  overcomplicated  geometries  are  permitted. 
In  addition,  the  calculations  can  be  efficiently  performed  in  an  element -by-element  fashion  to  carry  out  the 
large-scale,  time-dependent,  nonlinear  compulations  on  supercomputers. 

In  the  specific  problem  examined  in  this  paper,  the  changes  in  the  geometry  of  the  problem  require  a 
formulation  that  can  handle  moving  boundaries  and  interfaces.  The  Deformable-Spatial-Domain/Stabilized- 
Space-Time  (DSD/SST)  formulation  introduced  by  the  team  at  the  University  of  Minnesota  (UM)  in  1990 
has  this  capability.  This  is  an  accurate,  general-purpose  stabilized  finite  element  formulation  for 
computation  of  unsteady  flow'  problems  involving  free  surfaces,  two-liquid  interfaces,  and  fluid-structure 
and  fluid-particle  interactions.  In  the  DSD/SST  method,  the  stabilized  finite  element  formulation  of  the 
governing  equations  is  written  over  the  space-time  domain  of  the  problem,  and  therefore  the  deformation  of 
the  spatial  domain  with  respect  to  time  is  taken  into  account  automatically.  With  the  advanced 
stabilization  techniques  used  in  the  DSD/SST  formulation,  the  numerical  stability  challenges  involved  in 
the  problem  are  overcome  with  minimal  numerical  dissipation  and,  therefore,  with  minimal  loss  of 
accuracy. 


An  arbitarv  Lagrangian-Eulcrian  finite  element  method 
for  fluid-structure  interaction  problem 


Akira  ANJU,  Akira  MARUOKA  and  Mutsuto  KAWAHARA 
Department  of  Civil  Engineering,  Chuo  University,  Tokyo,  Japan 
1.  INTRODUCTION 

This  paper  presents  a  finite  element  analysis  of  a  fluid-structure  interaction  problem,  in  which  the  fluid  is 
treated  as  incompressible  viscous  flow,  and  a  structure  is  idealized  as  a  rigid  body  supported  by  elastic 
springs.  The  arbitrary  Lagrangian-Eulerian  (ALE)  method[t,2]  is  employed  to  solve  the  flow  field  around 
the  stmeture,  and  the  fractional  step  method[2,3]  is  adopted  for  the  time  integration.  For  the  numerical 
example,  the  present  method  is  applied  to  the  flow  analysis  around  the  oscillating  rectangular  cylinder. 
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1.  INTRODUCTION 

Massively  parallel  finite  element  simulations  of  incompressible  flows  with  the  Connection  Machine 
CM-5  and  CM-200  are  presented.  These  simulations  include  flow  problems  involving  moving  boundaries 
and  interfaces-such  as  fluid-body  and  fluidstructure  interactions— and  3D  flows.  For  problems  involving 
moving  boundaries  and  interfaces  we  use  the  deformable-spatial-domain/stabilized-space-time  formulation 
which  was  first  introduced  by  Tezduyar  et  al  (  1  ],  and  which  automatically  takes  into  account  the 
deformation  of  the  spatial  domain  with  respect  to  time.  To  update  the  mesh  as  the  spatial  domain  deforms 
with  time  in  this  type  of  problems,  special  mesh  moving  techniques  as  well  as  automatic  mesh  moving 
methods  are  used  to  minimize,  and  in  some  cases  totally  eliminate,  the  need  for  remeshing  (i.e.,  the  need 
for  generating  a  new  set  of  nodes  and  elements).  Remeshing  involves  projection  errors,  and  in  3D  problems 
the  cost  of  repetitive  mesh  generation  and  parallelization  set  up  (following  each  remeshing)  could  be 
prohibitive.  For  these  reasons,  we  find  mesh  update  strategies  which  minimize  the  frequency  of  remeshing 
desirable.  If  the  problem  does  not  involve  any  moving  boundaries  or  interfaces,  then,  to  avoid  the  cost 
associated  with  space-time  formulations,  we  prefer  to  use  the  pressure-stabilizing/Petrov-Galerkin  (PSPG) 
formulation  [2],  which  can  be  cast  in  the  context  of  a  semi-discrete  formulation.  The  stabilized  nature  of 
our  formulations  allows  us  to  use  equal-order  interpolation  functions  for  velocity  and  pressure. 

We  started  implementing  our  incompressible  flow  formulations  on  the  Connection  Machines  in  the 
fourth  quarter  of  1991.  We  reported  our  first  results  in  Tezduyar  et  al.  [3],  and  we  were  soon  able  to  solve  a 
large  class  of  practical  flow  problems,  such  as  fluid-body  and  simple  fluid-structure  interaction  problems 
[4]  and  3D  problems  [5-6],  In  all  these  computations  the  parallel  implementation  is  based  on  the 
assumption  that  the  mesh  is  unstructured,  and  this  gives  us  a  parallel  performance  with  very  little 
dependence  on  the  complexity  of  the  geometry'  and  mesh.  Furthermore,  all  these  computations  are  based  on 
implicit  time-integration  methods. 

These  new  computational  capabilities  are  now  allowing  us  to  solve  problems  which  we  were  not  able 
to  attempt  before.  Today,  essentially  all  our  computations  are  carried  out  on  the  massively  parallel 
platforms  of  the  Connection  Machines,  with  computational  speeds  which  are  two  orders  of  magnitude 
higher  than  what  was  available  to  us  two-three  years  ago  on  the  traditional  vector  supercomputers.  Today 
we  routinely  solve  timedependent,  3D  incompressible  flow  problems  and  even  3D  problems  with  moving 
boundaries  and  interfaces. 
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1.  INTRODUCTION 

Finite  element  computations  of  shallow  water  flows  and  contaminant  transport  can  be  applied  to  many 
practical  problems:  design  of  river,  coastal  and  offshore  structures,  disaster  prediction  and  other 
applications  related  to  hydrodynamic,  thermal  and  chemical  transport  behavior  in  oceans,  lakes,  and  rivers. 
In  this  context,  the  finite  element  method  is  applicable  to  complicated  water  and  land  configurations.  In 
practical  computation  of  this  type  of  problems,  it  is  essential  to  use  methods  which  are  as  efficient  and  fast 
as  the  available  hardware  allows.  Also,  in  this  type  of  problems,  computations  need  to  be  carried  out  over 
long  time  durations  to  properly  simulate  and  predict  the  phenomena  of  interest. 

In  recent  years,  massively  parallel  finite  element  computations  have  been  successfully  applied  to 
several  large-scale  compressible  and  incompressible  flow  problems,  including  those  involving  moving 
boundaries  and  interfaces  and  those  in  3D  [  1].  These  computations  demonstrated  the  availability  of  a  new 
level  of  finite  element  computational  capability  to  solve  practical  flow  problems.  With  the  need  for  a  high- 
performance  computing  environment  to  cany'  out  simulations  for  practical  problems  in  shallow  water  flows 
and  contaminant  transport,  in  this  paper  we  present  and  employ  a  parallel  explicit  finite  element  method  for 
computations  based  on  unstructured  grids.  The  finite  element  discretizations  are  based  on  a  three-step 
explicit  formulation  both  for  the  shallow  water  equations  and  the  advection-diffusion  equation  governing 
the  contaminant  transport.  In  these  discretizations,  for  numerical  stabilization,  we  use  selective  lumping  [2] 
for  the  shallow  water  equations  and  the  streamline-upwind/Petrov-Galerkin  (SUPG)  technique  [3]  for  the 
advectiondiffusion  equation.  Parallel  implementation  of  these  unstructured  grid-based  formulations 

carried  out  on  the  Connection  Machine  CM-5.  As  an  example,  we  carry  out  simulation  of  the  effect  of 
tidal  waves  on  the  Tokyo  Bay  and  the  spread  of  a  pollutant  injected  into  the  Tokyo  Bay. 
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1.  INTRODUCTION 

A  parallel  implementation  of  an  implicit  finite  element  formulation  for  incompressible  fluids  on  a 
distributed-memory  massively  parallel  computer  is  presented.  Data  stmctures  as  well  as  data  decomposition 
and  data  communication  algorithms  designed  for  distributed-memory  computers  are  discussed  in  the 
context  of  high  level  language  constructs  from  High  Performance  Fortran.  The  discussion  relies  primarily 
on  abstract  features  of  the  hardware  and  software  environment  and  should  be  applicable,  in  principle,  to  a 
variety  of  distributed-memory  systems.  The  implementation  is  carried  out  on  a  Connection  Machine  CM-5 
system  with  high  performance  communication  functions. 
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1.  INTRODUCTION 

Sloshing  is  a  free  surface  flow  problem  which  occurs  in  a  tank  subjected  to  forced  oscillation.  The 
main  difficulty  of  computation  in  large  amplitude  sloshing  analysis  arises  from  the  fact  that  both  the 
calculations  of  an  unknown  free  surface  position  and  the  free  surface  condition  should  be  satisfied 
simultaneously,  and  in  addition  the  fluid  is  strongly  distorted.  To  overcome  this  point,  a  new  Arbitrary 
Lagrangian-Eulerian  finite  element  method  for  three-dimensional  free  surface  flow  problems  in  a  tank  was 
developed. 
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Abstract 


Stabilized  finite  element  methods  such  as  SUPG  and  Galerkin/ieast  squares  (GLS) 
have  proved  to  be  effective  technologies  for  solving  the  Navier-Stokes  equations  and  general 
advective-diffusive  systems.  The  methods  were  initially  developed  for  incompressible  flows  in 


[1]  and  later  extended  to  compressible  flows  [25],  Our  work  in  compressible  flows  has 
emphasized  the  use  of  entropy  variables.  Although  the  method  based  on  entropy  variables 
possesses  unique  properties,  it  is  very  appealing  to  extend  the  formulation  to  other  sets  of 
variables,  such  as  primitive  variables.  Hughes  and  Tezduyar  [6],  Hansbo  and  Johnson  [7] 
and  Aliabadi,  Ray  and  Tezduyar  [8]  have  already  developed  similar  methods  within  the 
context  of  conservation  variables.  Possible  benefits  of  the  availability  of  the  various 
formulations  include  easier  implementation  of  the  GLS  finite  element  method  within  existing 
codes  and  extension  of  the  range  of  turbulence  models  that  can  be  used.  The  use  of 
primitive  variables  or  entropy  variables  also  allows  the  solution  of  incompressible  and 
compressible  flows  within  the  same  program. 

The  success  of  stabilized  finite  element  methods  for  the  laminar  Navier-Stokes 
equations  has  motivated  extension  to  turbulent  flows.  We  consider  Reynoldsaveraged 
systems  which  include  a  turbulent  kinetic  energy  transport  equation  ("oneequation  models"). 
We  identify  a  turbulent  entropy  production  (in)equality  associated  with  such  systems  and 
show  that  it  may  be  automatically  embedded  within  our  numerical  formulation,  thus 
incorporating  the  fundamental  nonlinear  stability  condition  ab  initio. 

We  also  present  a  few  numerical  calculations  illustrating  the  preceding  ideas.  In 
addition,  we  present  applications  on  massively  parallel  computers.  This  area  is  one  of  the 
most  exciting  in  the  history  of  scientific  computation.  The  engineering  potential  for  such 
machines  is  truly  mind  boggling.  We  briefly  illustrate  issues  of  load  balancing  via  a  parallelized 
domain  decomposition  algorithm,  and  present  timings  for  large  meshes  on  a  contemporary 
platform,  the  Connection  Machine  CM-5. 
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1  Introduction 

The  most  important  problem  of  computational  fluid  dynamics  is  the  accurate  prediction  of  turbulent 
phenomena  of  advection-dominated  flows  at  high  Reynolds  numbers.  In  order  to  accomplish  this 
prediction,  high  resolution  schemes  have  been  developed  in  each  field  of  finite  difference,  finite  volume 
and  finite  element  methods,  for  examples,  the  third-order  upwind  method  and  QUICK  by  Leonard  (1979) 
in  the  finite  difference  field,  MUSCL  by  Van  Leer  (1979)  and  SIMPLE  by  Patanker  (1980)  in  the 
volume  field,  SUPG  by  Hughes  et  al  (1977),  BTD  by  Dukowicz-Ramshow  (1979),  GSMAC  by  Tanahashi 
et  al.(1990),  upstream-downstream  selective  elements  method  by  Tabata-Fujita  (1991)  and  exponential-type 
Petrov-Galerkin  finite  element  method  by  Kakuda-Tosaka  (1992).  However,  upwind  technology  in  the 
finite  element  field  compared  with  the  other  numerical  fields  has  not  been  advanced.  The  finite  analytic 
method  by  Chen  (1987)  in  the  finite  difference  field  has  gotten  attractive.  In  this  method,  local  exact 
solution  of  steady  advection-  diffusion  equation  is  used  in  order  to  simplify  algorithms  and  to  obtain  the 
stability  at  high  Reynolds  numbers. 

The  purpose  of  this  paper  is  to  extend  the  finite  analytic  method  in  the  finite  difference  field  to  the 
finite  element  field.  This  possibility  is  attained  to  use  the  dual  spaces  composed  of  adjoint  operators. 
Accurate  interpolation  function  and  upwind-weighting  function  are  made  at  the  same  time  in  each  dual 
spaces.  Here  we  introduce  a  discrete  del  operator  and  an  adjoint  discrete  del  operator  in  the  dual  spaces. 
Those  operators  reduce  the  capacity  of  diffusion  matrix  and  get  the  high  speed  calculation  in  the  finite 
element  method.  The  basic  idea  of  stability  is  the  element  average  of  continuity  equation  and  the  node 
average  of  momentum  equation.  This  idea  in  the  unstructured  mesh  is  a  generalization  of  the  staggered 
mesh  in  the  MAC  method.  That  is  to  say,  the  conservation  of  mass  and  momentum  is  done  in  the  different 
domain.  After  mixing  the  two  equations  in  the  difference  domain,  we  will  get  the  numerical  solution  of  the 
incompressible  Navier-  Stokes  equation.  In  this  paper,  a  fundamental  theory  of  discretization  for  finite 
element  methods  will  be  discussed. 
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Third-Order  Upwind  Finite  Element  Simulation  of  Flow  around  Twin 

Square  Cylinders 
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1.  INTRODUCTION 

In  this  paper  numerical  results  of  flow  around  twin  square  cylinders  are  presented  on  the  basis  of 
numerical  simulation  of  the  incompressible  Navier-Stokes  equations  using  a  third-order  upwind  finite 
element  scheme[l].  It  is  well  known  that  a  flow  pattern  around  the  twin  square  cylinders  is  more 
complicated  than  a  flow  pattern  around  a  square  cylinder  because  of  interference  between  the  Karman 
vortices  behind  the  twin  square  cylinders.  We  therefore  show  profiles  of  the  interference  of  the  Karman 
vortices  occurred  behind  the  twin  square  cydinders. 
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element  scheme  based  on  the  choice  of 
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1.  INTRODUCTION 

In  [1]  we  have  developed  an  upwind  finite  element  scheme  for  the  incompressible  Navier-Stokes 
equations  at  high  Reynolds  numbers.  The  idea  of  the  upwind  technique  is  based  on  the  choice  of  up-  and 
down-wind  points  (see  Fig.l).  We  refer  to  [2,3]  for  the  details  and  the  advantages  of  this  scheme.  Here  we 
only  state  that  the  scheme  has  a  potential  to  approximate  the  convection  term  in  third-order  accuracy.  Thus, 
it  can  capture  the  behavior  of  high  Reynolds  number  flows  well. 

In  this  paper  we  present  some  techniques  for  high  performance  implementation  of  the  scheme  in  terms 
of  computation  speed  and  memory  size. 
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1.  INTRODUCTION 

We  present  a  brief  description  of  a  finite  element  Solver  implemented  on  the  Cormection  Machine 
CM-5  system.  A  more  detailed  presentation  of  the  issues  in  such  an  implementation  can  be  found  in  [1,  2], 
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1.  INTRODUCTION 

In  this  paper,  we  present  a  finite  element  solution  procedure  for  solving  large  scale  time- 
dependent  incompressible  viscous  flow  problems.  It  is  a  velocity-explicit/pressure-implicit  time¬ 
marching  algorithm.  For  memory-saving,  the  Qj-Po  elements  are  employed  and  the  discrete 
pressure  Poisson  equation  is  consistently  derived  from  the  mixed  finite  element  formulation. 
The  matrix  system  is  solved  by  an  element-byelement  conjugate  gradient  method,  which  is 
cost-effective  in  both  storage  and  CPU  time. 
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1.  INTRODUCTION 

Parallelization  of  a  flow  solver  using  a  machine  independent  data  base  program  GPAR  [  1], 
developed  primarily  for  parallelization  of  CFD  programs,  is  presented.  In  this  approach  the  solution  block 


domain  is  subdivided  into  a  number  of  solution  blocks  greater  than  or  equal  to  the  number  of  available 
computers.  The  blocks  are  distributed  to  available  computers  on  a  network  to  solve  the  same  problem  in 
parallel.  The  data  base  program  handles  blocks  and  their  interface  information.  The  machine  independence 
is  achieved  by  utilizing  a  machine  portable  parallel  library  APPL  [2]  developed  at  NASA  Lewis  Research 
Center.  The  flow  solver  consists  of  separate  block  and  interface  solvers  [3].  A  given  computational  grid  is 
subdivided  into  blocks  using  an  automatic  divider  algorithm  which  prepares  the  block  and  interface 
information  to  GPAR.  An  explicit  time  integration  scheme  is  used  following  an  upwinded  finite  element 
formulation  of  unsteady  Euler  equations. 

Since  for  accurate  calculation  of  lift  and  pressure  coefficients  around  aerodynamic  configurations 
large  size  problems  must  be  solved  for  many  time  steps,  memory  and  CPU  requirements  of  unsteady  flow 
solvers  are  rather  excessive.  The  parallel  approach  presented  herein  provides  a  practical  utilization  of 
networked  workstations  for  this  puqwse. 
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1.  INTRODUCTION 

The  cellular  automaton  method  (the  CA  method)  is  a  particle  simulation  method  for  an  incompressible 
fluid  [1][2].  Recently,  some  complicated  flows  with  multi  components  have  been  simulated  through  the 
microscopic  approach  of  the  CA  method  (e.g.  immiscible  two  phase  flow  [3],  share  layer  analysis  [4]). 

On  the  other  hand,  it  is  known  that  the  discrete  structure  of  the  CA  method  is  suited  for  parallel 
processing.  In  fact,  we  have  developed  the  parallel  CA  method  with  1-dimensional  domain  decomposition 
technique  and  have  obtained  the  speedup  factor  of  more  than  63  by  a  Fujitsu  APIOOO  with  64  processors 
in  327680  node  calculation  [5].  The  high  efficiency  of  the  parallel  CA  method  will  be  maintained  in 
massively  parallel  processing  (MPP)  for  the  practical  use. 

In  this  paper,  Poiseuille  flow  and  immiscible  two  phase  flow  simulation  by  the  parallel  CA  method 
are  presented  to  estimate  the  suitability'  for  the  MPP. 
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The  analyses  of  buoyancy -driven  flows  in  fluid  layers  heated  from  below  has  been  the  of  many  theoretical 
and  experimental  investigations  since  the  observations  of  Tyndall  in  1863  [1]  and  the  analysis  of  Lord 
Rayleigh  in  1916  [2],  Such  fundamental  aspects  as  the  onset  of  convective  instability,  transition  to 
periodic  and  non-periodic  motions,  influence  of  the  form  of  the  perturbation  and  boundary  conditions,  etc. 
have  been  extensively  considered  via  linear  and  weakly  nonlinear  analyses.  More  recently  attention  has 
again  been  focused  on  the  effect  of  rigid  lateral  boundaries  on  the  threshold  of  stability  and  the  form  of  the 
ensuing  flow;  in  particular,  rectangular  and  circular  cylindrical  domains  have  been  investigated  over  a  range 
of  aspect  ratio.  Aside  from  the  obvious  fundamental  aspects  of  such  studies,  insight  is  also  gleaned  into 
possible  convection  in  the  vapor,  or  liquid,  nutrient  phase  in  certain  crystal-growing  techniques  as,  for 
example,  the  Bridgman  technique  (Carruthers  (3]).  The  interest  here  is  primarily  in  circular  cylindrical 
domains  in  both  a  1-g  terrestrial  as  well  as  a  low-g  near-earth  orbit  where  the  residual  gravity  vector  can 
vary  both  in  magnitude  and  direction  as,  for  example,  in  the  g-jitter  phenomena.  (See,  for  example,  Ostrach 

[4]). 

A  number  of  studies  devoted  to  convection  in  vertical  confined  cylinders  heated  from  below  have 
appeared  in  the  literature  but  only  a  few  have  dealt  with  the  rigid  (no-slip)  wall  case  and  only  a  very 
limited  number  of  fully  three-dimensional  numerical  simulations  are  available.  (See,  for  example,  Smutek 
et  al.  [5]  for  gases,  Pr=0.7  and  Neumann  [6]  for  liquids,  Pr=6.7).  The  previous  theoretical  studies  were 
limited  to  linear  and  weakly  nonlinear  stability  analyses.  (See  Hardin  et  ai.  [7]  and  Hardin  and  Sani  [8]  for 
a  summary  and  new  results.)  While  these  analyses  establish  that  the  onset  and  slightly  supercritical  state  of 
steady  convection  is  strongly  dependent  on  aspect  ratio  (Charlson  and  Sani  [9]  and  Buell  and  Cation  [10]), 
they  are  silent  about  the  flow  and  structure  at  elevated  supercritical  states.  The  current  study  is  focused  on 
using  transient,  3D  numerical  simulations  to  confirm  previous  results  as  well  as  to  explore  the  supercritical 
flow  and  structure  in  such  sj'slems. 
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INTRODUCTION 

For  effective  control  of  the  indoor  climate  in  an  atrium,  it  is  very  important  to  predict  the  and  flow 
fields  during  design.  One  of  the  methods  for  this  prediction  is  numerical  simulation.  We  have  developed 
the  numerical  simulation  system,  which  consists  of  three  subsystems  as  follows;  (1)  simulation  system,  of 
the  direct  solar  radiation,  tJie  diffused  solar  radiation  and  the  diffused  reflection  of  the  solar  radiation,  (2) 
simulation  system  of  the  radiative  heat  transfer,  (3)  simulation  system  of  the  temperature  and  flow 
distribution  based  on  k  -  e  model  of  turbulence  solved  by  GSMAC3D-FEM[1,2]. 

On  the  other  hand,  to  verify  the  accuracy  and  the  reliability  of  those  systems,  comparison  between 
results  from  expieriments  and  numerical  analysis  was  carried  out  in  the  test  atrium  for  experiments.  As  a 
result  of  comparison,  it  became  obvious  that  results  of  numerical  analysis  agreed  with  those  of  experiments 
by  coupling  the  calculation  of  radiative  heat  transfer  and  that  of  convective  heat  transport. 
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1.  INTRODUCTION 

State-of-the-art,  high-performance  computing  has  proven  to  be  an  invaluable  tool  for  addressing  a 
spectmm  of  issues  in  mani^acturing  technology.  Understanding  the  complex  physical  phenomena  involved 
with  manufacturing  processes  poses  challenges  comparable  to  the  "Grand  Challenge"  supercomputer 
problems  of  the  pure  physical  sciences  and  is,  perhaps,  even  more  compelling  due  to  the  obvious  economic 
import  of  manufacturing  activity. 

Our  interests  center  on  processes  employed  to  produce  advanced  materials,  such  as  large,  single 
crystals  for  photonic  and  electronic  devices  [1  3]  and  ceramic  materials  [4,5].  Such  processes  are 
characterized  by  complex  geometries  which  are  often  time-dependent,  coupled  incompressible  flows,  heat 
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and  mass  transfer,  chemical  reaction,  and  phase  change.  As  the  technological  complexity  of  such  processes 
continues  to  grow,  increasingly  powerful  techniques  must  be  employed  for  solution  of  the  mathematical 
models  which  describe  their  physical  behavior  [6], 

Here,  we  report  on  continuing  progress  in  developing  finite  element  methods  based  on  algorithms 
which  perform  effectively  on  distributed  memory,  massively  parallel  supercomputers,  such  as  the  Thinking 
Machine  Corporation  CM-5.  Our  initial  efforts  were  reported  in  [7], 
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1.  Introduction 

A  numerical  method  is  presented  to  analyze  a  transient  convective  diffusion  problem  with  a  first-order 
chemical  reaction  defined  on  an  infinite  domain.  The  present  method  for  transient  problems  is  based  on  the 
combined  boundary  and  finite  element  method[l]  using  Laplace  Transform.  For  the  time  integration,  an 
inverse  Laplace  transform  presented  by  Hosono[2]  has  been  successfully  employed. 


The  finite  element  method  is  one  of  the  useful  numerical  tools  to  analyze  these  problems.  It  is 
however,  difficult  to  analyze  a  problem  defined  on  the  infinite  domain  by  use  of  the  finite  element 
methods.  The  boundary  element  method  satisfies  unconditionally  the  infinite  boundary  condition  by 
assuming  that  the  problem  is  linear.  However,  the  boundary  element  method  has  difficulties  in  solving 
problems  involving  non-homogeneous  fields  because  it  is  difficult  to  obtain  the  fundamental  solution 
explicitly.  If  we  appropriately  combine  the  finite  element  and  the  boundary  element  methods,  it  is  obvious 
that  the  method  which  covers  the  difficulties  of  both  methods  can  be  obtained.  For  steady  problems  the 
combined  method  is  recognized  as  successful  numerical  technique  to  deal  with  the  infinite  domain.  This 
method  is  also  useful  in  transient  problems.  In  the  computation  of  the  transient  problem  using  the 
boundary  element  method,  the  domain  discretization  should  be  employed  to  integrate  over  the  infinite 
domain.  Furthermore,  difficulties  arise  when  the  integral  domain  should  be  expanded  according  to  the 
expansion  of  the  diffusion  domain.  The  above-mentioned  difficulties  can  be  overcome  by  introducing  the 
Laplace  transform.  The  time  dependence  of  the  problem  can  temporarily  be  removed  by  this  method. 
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1.  INTRODUCTION 

Computational  mechanics  technolog)'  has  progressed  dramatically.  The  use  of  high  performance 
computers  such  as  massively  parallel  computers  has  enabled  us  to  analyze  nonlinear  large  scale  problems  in 
various  engineering  fields.  Generally,  most  of  the  nonlinear  analyses  involve  several  parameters,  which 
greatly  influence  both  the  numerical  results  and  the  stability.  These  parameters,  e  g.,  the  penalty  parameter 
for  incompressible  viscous  flow  analysis  and  the  time  step  of  transient  analysis  are  to  be  determined 
appropriately  and  dynamically  referring  to  convergence  behaviors.  However,  since  nonlinear  analysis  in 
each  field  has  its  own  empirical  knowledge  for  choosing  optimal  parameters,  it  is  difficult  to  decide  the 
general  optimization  rule  or  the  procedure  for  nonlinear  finite  element  analyses. 

The  present  paper  describes  a  new  methodology  for  optimizing  multiple  parameters  of  nonlinear  finite 
element  analyses  using  the  hierarchical  neural  network.  To  demonstrate  the  performance  of  the  present 
procedure,  it  is  applied  to  the  parameter  optimization  of  the  augmented  Lagrangian  method  for  steady  state 
analysis  of  incompressible  viscous  flow  and  the  time  step  optimization  of  pseudo  time-dependent  stress 
analysis  for  incompressible  inelastic  stmeture. 
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1.  INTRODUCTION 

Computer  codes  for  wave  propagation  and  impact  have  matured  considerably  since  their  initial 
development  some  35  years  ago.  They  now  sen'e  as  valuable  tools  in  studies  of  materials  and  structures 
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subjected  to  intense  impulsive  loading.  Numerical  simulations  of  high-velocity  impact  phenomena  in  two 
dimensions  have  been  performed  routinely  for  a  number  of  years.  The  advent  of  high-performance  vector 
and  parallel  computers  has  made  three-dimensional  solutions  of  the  governing  equations  tractable. 
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1.  INTRODUCTION 

In  this  paper,  we  shall  focus  on  the  efficient  use  of  modem  computing  equipment  for  high-perfomiance 
finite  element  computations.  We  devote  our  attention  to  the  evaluation  of  computational  strategies,  data 
stmcture  and  flow,  algorithmic  design  and  programming  techniques.  In  particular,  we  investigate  the 
effective  use  of  networked  workstations  equipped  with  high-performance  microprocessors.  Such  systems 
often  outperform  multi-user  supercomputing  facilities  in  terms  of  the  elapsed  time  to  complete  certain 
computationally  intensive  tasks.  Our  discussion  will  be  restricted  to  the  implicit  numerical  processes 
involving  the  formulation  and  factorization  of  finite  element  coefficient  arrays. 

The  finite  element  method  has  traditionally  been  deployed  on  the  most  powerful  equipment 
because  of  its  computational  intensity.  Progress  in  hardware  and  software  technologies  has  made  it  possible 
to  tackle  complicated  engineering  and  scientific  problems  on  powerful  desktop  systenis.  Problems  only 
tractable  on  supercomputers  in  the  recent  past  are  now  routinely  solved  using  these  systems. 

The  use  of  massively  parallel  systems  are  presently  limited  to  the  research  community.  Many 
architecture  and  programming  paradigms  are  being  proposed  and  evaluated.  Considerable  time  will  be 
required  for  such  systems  to  be  commonly  accepted  by  engineering  community  and  used  for  practical 
applications.  Coarse  grain  parallelization  such  as  symmetric  multi-processor  systems  are  emerging  on 
workstations  and  servers.  We  would  focus  our  future  investigations  on  the  effective  use  of  such  systems  as 
a  natural  extension  to  what  is  presented  in  this  paper. 
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1.  INTRODUCTION 

The  conventional  finite  element  method  based  on  the  piecewise  polynomial  interpolations  has  limited 
capability  to  analyze  the  phenomena  in  the  wholeor  half-space.  Here  infinite-domain  interpolations  which 
can  be  conformably  connected  with  the  isoparametric  elements  are  of  great  significance. 

In  such  infinite  element,  our  trial  function  is  to  be  selected  over  each  (semi-)infinite  domain  so  that 
the  dominant  decaying  mode  is  reproducible,  at  least  approximately,  in  the  global  system.  Some  infinite 
elements  adopt  the  special  non-polynomial  trial  function  space  under  the  usual  polynomial  mapping  [13, 
while  the  others  depend  on  the  specific  mapping  called  the  infinity  mapping  with  polynomial 
interpolations  [23,  [3], 

This  paper  is  devoted  to  establish  the  unified  theory  in  the  concerning  infinite  elementology. 
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Maebashi  City  College  of  Technology,  Gumma,  JAPAN 
1  INTRODUCTION 

This  paper  presents  a  computational  technique  of  the  moving  water  domain  based  on  shallow  water 
equations  by  two-step  explicit  finite  element  method.  It  is  impiortant  to  know  the  speed  of  water  up  on  the 
drying  bed,  when  the  dam  is  broken.  The  water  boundary  is  calculated  by  moving  boundary  techmque 
based  on  finite  element  method.  In  this  study,  the  moving  boundary  technique  is  used  by  finding  if  an 
element  has  water  depth  or  not. 

It  would  be  better  to  use  the  same  technique  to  find  the  broken  dam  problem  as  in  the  flood  flow  river. 
However,  some  declinations  occurred  using  such  technique,  it  is  convenient  to  operate  with  such  method. 
These  declinations  are  present  from  the  simulation  of  the  highest  jxiints  of  the  water  after  the  dam  is 
broken,  as  if  these  points  are  part  of  the  flood  flow  river.  In  this  case,  the  velocity  of  the  boundary  nodes  is 
zero.  But  in  the  broken  dam  problem,  the  velocity  of  these  points  is  not  equal  zero. 

In  this  paper,  this  problem  is  treated  and  new  technique  is  introduced  to  reduce  such  abnormal 
phenomenon.  The  numerical  results  of  a  simple  broken  dam  example,  are  found  by  old  and  new'  moving 
boundary  technique. 


Simulations  of  Coupled  Viscous  and  Porous 
Flow  Problems  ^ 

D.K.  Gartling,  C.E.  Hickox  and  R.C.  Givler 

Computational  Fluid  Dynamics  Department  Sandia  National  Laboratories 
Albuquerque,  New  Mexico  87185 


1.  INTRODUCTION 

Conjugate  problems  involving  a  viscous  flow  and  an  adjoining  fluid  saturated  porous  medium 
represent  an  important  class  of  problems  for  a  variety  of  technologies.  In  some  applications  the  interface 
between  the  porous  material  and  the  clear  fluid  is  distinct  and  its  location  is  known  a  priori.  Typical  of 
these  applications  are  packed  beds  for  thermal  or  chemical  processing,  insulating  structures,  filters  and 
porous  bearings.  A  more  difficult  problem  involves  a  free  boundary  interface,  where  the  location  of  the 
clear  fluid/porous  boundary  is  determined  by  the  variation  of  one  or  more  dependent  variables  and  therefore 
evolves  with  the  dynamics  of  the  flow.  The  prototypical  application  for  this  type  of  conjugate  problem  is 
found  in  the  solidification  of  metal  alloys  where  the  evolving  dendrite  structure  is  represented  by  a  porous 
medium.  The  most  general  situation  for  all  applications  is  one  in  which  the  interface  between  the  bulk 
fluid  and  the  porous  region  is  completely  open,  allowing  the  free  transport  of  fluid,  heat,  etc.  between  the 
regions. 

The  formulation  of  this  type  of  conjugate  flow  problem  requires  that  two  inter-related  questions  be 
resolved.  The  equations  describing  the  flow  of  the  viscous  fluid,  the  NavierStokes  equations,  are  well- 
established  and  are  not  debated  for  this  application.  The  appropriate  description  for  flow  in  the  porous 
medium  is,  however,  not  as  evident,  nor  is  the  related  issue  of  appropriate  boundary  and  interface 
conditions  clearly  understood.  Though  Darcy's  law  is  a  well  known  and  justifiable  description  for  saturated 
porous  flows,  it  has  been  criticized  for  its  range  of  applicability  and  the  limited  types  of  boundary  and 
interface  conditions  that  may  be  reliably  employed  with  the  model.  Non-Darcy  descriptions,  such  as  the 
Brinkman  equation,  seek  to  remedy  the  boundary  and  interface  condition  problems  through  the  addition  of 
higher  order  terms  to  the  flow  equation.  Though  this  approach  eases  the  boundary  condition  problem, 
additional  material  parameters  are  introduced  that  are  difficult  to  quantify.  The  problem  is  therefore  to 
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select  an  appropriate  equation  for  the  porous  flow  region  in  conjunction  with  the  proper  interface  conditions 
to  couple  it  to  the  bulk  flow.  An  additional  concern  involves  the  appropriate  method  of  transition  from  a 
viscous  to  a  porous  flow  when  the  interface  is  not  distinct  and  is  described  by  a  porosity  variation. 

The  objective  of  the  present  work  was  to  develop  a  general  method  for  solving  conjugate  fluid  and 
porous  flow  problems  using  a  finite  element  method.  In  pursuing  this  goal  several  combinations  of  porous 
flow  equations  and  interface  conditions  were  examined,  compared  and  evaluated.  The  following  section 
presents  the  general  mathematical  description  of  the  equations  used  in  the  study;  a  following  section 
outlines  a  number  of  non-Darcy  effects  and  possible  interface  conditions  that  have  been  evaluated. 

^This  work  was  performed  at  Sandia  National  Laboratories  and  is  supported  by  the  U.S.  Department 
of  Energy  under  contract  DE-AC04-94-AL85000. 
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1.  INTRODUCTION 

The  pressure-strain  correlation  term  in  the  differential  stress  model  is  the  most  important  term  for 
accurate  estimation  of  wind  load  on  buildings[ll.  The  larger  the  turbulent  kinetic  energy  on  the  inflow 
boundary  becomes,  the  more  significant  the  error  caused  by  the  imperfection  of  this  term  becomes. 
Moreover  the  cost-efficiency  of  the  differential  stress  model  is  inferior  to  the  k-e  turbulence  model. 
Therefore  most  of  the  practical  applications  of  turbulence  model  in  wind  engineering  field  are  executed  by 
the  k-e  turbulence  model. 

This  paper  presents  an  improved  k-e  turbulence  model  for  accurate  estimation  of  wind  load  on  bluff 
body.  The  deceleration  effect  through  impingement  against  wall  and  adverse  pressure  gradient,  and  the 
acceleration  effect  through  favorable  pressure  gradient  are  introduced  to  the  turbulent  kinetic  energy 
production  term  of  the  conventional  standard  k-e  turbulence  model.  The  flows  around  square  prism  in 
turbulent  streams  are  computed  for  validation  of  the  improved  model,  and  the  pressure  distribution  on  tlie 
surface  is  improved  dramatically  compared  witli  ilie  standard  k-e  turbulence  model. 
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A  numerical  study  on  vortex  shedding  around  a  heated  / 
cooled  circular  cylinder  by  three-step  Taylor-Galerkin 
method 

Katsumori  HATANAKA 
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1.  INTRODUCTION 

The  vortex  shedding  around  circular  cylinder  has  attracted  much  attention  from  researchers  who  work 
on  both  experimental  and  numerical  fluid  dynamics.  The  vortex  shedding  in  the  mixed  convection 
phenomeivon  is  more  complex  due  to  the  temperature  induced  buoyancy  effect  and  is  much  more  interesting 
than  an  ordinary  vortex  shedding  problem.  (See  [  1  ]  or  [2],  for  example.) 

In  this  paper  the  vortex  shedding  around  a  heated  /  cooled  circular  cylinder  has  been  investigated 
numerically  by  the  three-step  Taylor-Galerkin  /  velocity  correction  method  [3].  The  attention  of  the  present 
paper  is  directed  to  the  change  of  the  vortex  shedding  patterns  under  influence  of  buoyancy  forces  at  a 
constant  Reynolds  number.  The  flow  problem  under  consideration  has  been  successfully  simulated  by  the 
present  finite  element  scheme. 
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Applications  of  CFD  to  Wind  Engineering 
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l.introduction 

We  are  preparing  four  topics  based  on  our  recent  researches  for  the  symposium.  Those  are  the  influence 
of  open  boundary  width  on  numerical  results  of  pressure  coefficients,  response  characteristics  of  a  vortex- 
exited  oscillation  of  cylinder  in  low  Reynolds  numbers  flow,  response  characteristics  of  a  vortex-exited 
oscillation  of  cylinder  in  high  Reynolds  numbers  flow  and  a  term  wise  comparison  between  FDMs  and 
FEMs  by  using  Navier-Stokes  equatioa  In  this  paper,  a  brief  description  on  the  four  topics  is  described. 
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Analyses  Of  Mold  Flow  And  Residual  Stress  In  Resin  Transfer 
Molding  Using  The  Homogenization  Method 

Whie  Chang  and  Noboru  Kikuchi 

The  University  of  Michigan,  Ann  Arbor,  Michigan  48109 

Resin  Transfer  Molding 

Resin  transfer  Molding  (RTM)  is  a  relatively  new  process  for  manufacturing  fiber  reinforced 
polymer  composites.  In  RTM,  fiber  mats  which  act  as  reinforcements  are  cut  into  the  desired  shape  and 
arranged  inside  the  mold  cavity.  Then  the  mold  is  closed  and  the  resin  is  injected  through  the  injection 
ports  of  the  mold  progressively  impregnating  the  fiber  mats.  When  the  mold  is  completely  filled,  the  resin 
solidifies  in  the  curing  phase  in  which  chemical  reaction  and  significant  temperature  variation  take  place. 
After  the  curing  the  product  is  taken  out  of  the  mold. 

The  process  analysis  is  a  very  important  step  in  the  product  development  stages  of  RTM.  It 
consists  of  the  mold  flow  analysis  of  the  mold  filling  phase  and  the  residual  stress  analysis  in  the  curing 
phase.  The  former  can  provide  useful  information  to  a  mold  designer  such  as  pressure  distributions  and 
flow  front  profiles  within  the  mold  cavity  which  can  be  used  to  predict  operating  parameters,  possible 
location  of  defects,  hard  to  fill  regions,  and  proper  locations  of  air  tabs  and  injection  ports.  The  latter  can 
be  used  to  predict  the  residual  stresses  which  are,  in  some  cases,  large  enough  to  cause  micro  cracking  and 
exhaust  the  strength  of  the  material.  Conventionally,  these  data  were  obtained  based  on  experience  and 
repeated  experiments  which  are  time  consuming  and  inefficient.  In  order  to  optimize  the  manufacturing 
process  and  reduce  the  cost,  numerical  analysis  of  the  RTM  process  is  necessary. 

In  this  work,  a  complete  analyses  are  given  which  encompasses  the  formulations  and  the 
numerical  implementations  of  the  mold  flow  and  of  the  residual  stress.  The  mold  flow  analysis  deals  with 
the  resin  flow  through  the  complex  pore  structures  of  the  reinforcing  fiber  perform,  and  the  residual  stress 
analysis  deals  with  the  stress  distribution  in  the  micro  structure  of  the  composite.  In  both  analyses,  the 
homogenization  method  is  utilized  extensively  due  to  the  complexity  of  the  geometry. 

The  homogenization  theory  deals  with  the  partial  differential  equations  of  physics  in  heterogeneous 
materials  with  periodic  structures  when  the  characteristic  length  of  the  period  is  small.  In  the 
homogenization  theoiy',  it  is  assumed  that  the  composite  is  locally  formed  by  spatial  repetition  of  a 
"microscopic  cell".  Using  this  assumption,  an  equivalent  (or  "homogenized")  material  model  is  found  by 
solving  the  behavior  within  the  microscopic  cell.  The  theory  appeared  in  the  1970's  and  has  been  the 
subject  of  considerable  research  in  different  areas  of  applied  mechanics.  The  application  of  the  theory  to  the 
analyses  of  composite  materials  can  be  found  in  Sanchez-Palencia  [4]  and  some  computational  works  can 
be  found  in  Guedes[81  and  Cheng[3].  The  fundamentals  of  the  theory  can  be  found  in  the  works  of 
Lions[7],  Benssousan,  Liorrs,  and  Papanicoulau  [1],  and  01einik[9]. 
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1.  INTRODUCTION 

It  is  quite  useful  to  develop  an  adaptive  mesh  rezoning  technique  which  adjusts  to  the  moving  domain 
and  update  the  mesh,  yet,  maintaining  the  element  connectivity,  since  the  Space-time  Finite  Element 
Method  is  considered  as  a  promising  framework  for  fluid-structure  interaction  (FSI)  problems.  We  have 
developed  the  adaptive  mesh  rezoning  technique  with  the  least-square  constrained  conditions,  which  save 
the  element,  locating  close  to  the  moving  boundary',  from  breakdown  trouble  and  keep  well-conditioned 
shape.  The  our  technique  is  entirely  general  in  that  it  can  be  applied  to  structured  and  unstructured  meshes, 
effectively. 
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Adsorption,  in  which  molecular  constituents  of  a  fluid  are  partitioned  by  virtue  of  having  different 
affinities  for  a  solid,  is  used  to  separate  bulk  mixtures  as  well  as  remove  trace  impurities.  Adsorption 
processes  are  typically  carried  out  in  large  columns  packed  with  a  multitude  of  small  porous  pellets  to 
promote  contact  between  fluid  and  solid.  Species  in  the  feed  stream  are  fractionated  by  the  solid  on  an 
'adsorption'  step  which  is  then  followed  by  a  sequence  of  'rinse'  or  'purge'  steps  to  clean  out  the  solid  so 
that  it  may  be  used  again.  The  dynamics  vary'  in  detail  throughout  the  column,  but  at  each  point  eventually 
a  cyclic  or  periodic  state  is  reached  in  which  the  particles  undergo  equal  amounts  of  adsorption  and 
desorption. 
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1.  INTRODUCTION 

A  large  storage  requirement  of  the  finite  element  method,  which  is  a  result  of  using  unstructured  meshes 
and  the  need  to  solve  the  system  of  equations.  The  one-point  quadrature  technique  proposed  by  Gresho  et 
al.[l]  is  known  to  be  quite  efficient  in  saving  computational  storage.  By  assuming  the  one-point 
quadrature,  the  diffusion  and  the  advection  matrices  become  storage-free  since  they  can  be  computed  from 
the  gradient  Matrix.  Some  numerical  dumping  needs  to  be  introduced  to  suppress  physically  meaningless 
oscillations  called  the  'hourglass  mode',  which  might  be  caused  by  the  reduced  (  one-point )  integration  of 
the  diffusion  matrix.  The  authors  have  developed  an  hourglass  controller[2],  which  can  incorporate  the 
effect  of  an  aspect  ratio  of  element.  The  authors  also  have  proposed  a  matrix-storage  free  formulation.  Since 
the  gradient  matrix  still  occupies  large  amount  of  memory,  by  virtue  of  developing  an  efficient  algorithm 
to  compute  the  gradient  matrix,  a  matrix-storage  free  formulation  has  been  made  possible.  Comparably 
practical  problems  having  several  tens  of  thousands  of  nodes  can  be  easily  handled  by  currently  available 
engineering  workstations. 
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1.  INTRODUCTION 

Multiple  scale  methods  [1-4],  which  are  based  on  discrete  and  continuous  reproducing  kernels,  wavelets, 
and  integral  window  transforms,  are  developed.  In  this  development,  a  microscope  is  constructed  with  a 
flexible  space/time  localized  window  function  which  translates  and  dilates  in  space/time  to  cover  the  entire 
domain  of  interest.  This  microscope  can  magnify,  examine,  and  record  the  image  of  the  various 
scales/frequencies  of  the  response  locally  within  the  support  of  the  window  function.  The  degree  of 
magnification  will  depend  on  the  power  of  the  microscope,  a  flexible  space-scale/time-frequency  window 
function.  This  complete  characterization  of  the  unknown  respionse  is  performed  through  the  integral 
window  transform.  This  localization  process  can  be  achieved  by  ilating  the  flexible  multiple-scale  window 
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function.  A  parameter  defined  by  the  energy  error  ratio  is  introduced  to  pin  point  the  scale/frequency 
domain  of  interest,  and  to  prevent  aliasing  from  the  sampling  procedure.  The  zoom  in  and  zoom  out 
capability  of  the  window  function  is  especially  useful  in  examining  complex  flow  phenomena,  such  as 
flow  induced  vibration,  dynamic  stability  of  flow-stmcture  interaction,  turbulence  structures,  and  high 
frequency  structural  dynamic  response.  Numerical  examples,  which  include  the  Helmholtz  equation  and  the 
advection-diffusion  equation,  are  presented  and  compared  with  finite  element  results. 
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1.  INTRODUCTION 

In  general,  water  quality  control  is  one  of  the  most  important  social  problems  especially  closed  water 
region.  If  we  want  control  this  water  pollution  problem,  we  have  to  consider  the  inequality  constrains. 
Until  now,  we  have  been  studying  about  the  dynamic  optimization  without  constrain  problem  for  first 
step.  This  type  problem  can  be  solved  by  Conjugate  Gradient  (CG)  Method,  Dynamic  Programming  (DP) 
Method,  etc..  However  these  methods  caimot  be  solved  under  inequality  constrain.  Therefore  this  paper 
proposed  a  method  based  the  Sakawa-Shindo  method  to  solve  the  water  pollution  control  problem  with 
inequality  constraints  on  functions  of  control. 

In  this  study,  the  optimal  control  analysis  has  been  applied  to  water  pollution  problem.  The  main 
purpose  of  control  of  water  pollution  problem  is  to  reduce  the  contaminated  water,  i.e.  the  determination  of 
control  values  to  keep  the  water  quality  which  satisfy’  a  criterion  of  concentration  at  arbitary  point  in  the 
flow  field. 

In  the  present  optimal  control  analysis,  the  linear  shallow  water  equation  and  the  convection  diffusion 
equation  are  used  as  the  governing  equations  for  flow  field.  The  explicit  Euler  method  is  applied  to 
discretize  the  shallow  water  equation  and  the  convection  diffusion  equation.  The  simulations  of  simple 
model  and  applications  to  the  water  pollution  problems  have  been  carried  out.  The  results  of  the  method 
are  quite  resortable.  The  application  for  practical  problem  will  be  carried  out. 
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1.  A  THREE-FIELD  FORMULATION  OF  COUPLED  AEROELASTIC  PROBLEMS 

In  order  to  predict  the  dynamic  response  of  a  rigid  or  flexible  structure  in  a  fluid  flow,  the  equations  of 
motion  of  the  structure  and  the  fluid  must  be  solved  simultaneously.  The  most  difficult  part  of  handling 
numerically  the  fluid/structure  coupling  stems  from  the  fact  that  the  structural  equations  are  usually 
formulated  with  material  (Lagrangian)  co-ordinates,  while  the  fluid  equations  are  typically  written  using 
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spatial  (Eulerian)  co-ordinates.  Therefore,  a  straightforward  approach  to  the  solution  of  the  coupled 
fluid/structure  dynamic  equations  requires  moving  at  each  time  step  at  least  the  portions  of  the  fluid  grid 
that  are  close  to  the  moving  structure.  This  can  be  appropriate  for  small  displacements  of  the  stmcture  but 
may  lead  to  severe  grid  distortions  when  the  structure  undergoes  large  motion.  Several  different  approaches 
have  emerged  as  an  alternative  to  partial  re-gridding  in  transient  aeroelastic  computations,  among  which  we 
note  the  Arbitrary  Lagrangian  Eulerian  (ALE)  formulation  [1],  the  co-rotational  approach  [2,3],  and 
dynamic  meshes  [4].  All  of  these  approaches  treat  a  computational  aeroelastic  problem  as  a  two-field 
problem. 
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1.  INTRODUCTION 

In  the  coastal  management  for  the  coastal  defence  against  erosion,  storm  surge,  etc,  and  for  the  use  as 
hatbouts  and  resort  areas,  it  is  very  important  to  estimate  or  predict  the  currents  such  as  tidal  or  near  shore 
currents,  wave  propagation,  tidal  elevation  and  topography  change.  However,  for  the  reason  that  these 
phenomena  are  very  complicated  and  stochastic,  it  is  difficult  to  estimate  them  determinately.  On  the  other 
hand,  sometimes  they  are  observed  or  monitored  using  various  kinds  of  equipment.  Those  data  often  in¬ 
cludes  erroneous  disturbance.  In  this  paper,  the  Kalman  Filtering  tectmiqe  is  applied  to  the  wave 
propagation  analysis,  consequently  it  is  possible  to  take  into  account  the  mesurement  data  efficiently  to  get 
much  information  about  the  wave  movement. 
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1.  INTRODUCTION 

We  consider  finite  element  approximations  to  axisymmetric  flow  problems.  The  differential  operators 
in  these  problems  have  coefficients  of  negative  powers  of  the  distance  from  the  axis.  So  we  have  to  treat 
weighted  Sobolev  spaces  [1].  The  radius  and  the  axis  components  of  the  velocity  belong  to  different 
function  spaces.  So  we  have  to  deal  with  three  weighted  Sobolev  spaces  for  the  velocity  and  the  pressure. 

As  the  model  problem  we  consider  an  axisymmetric  stationary  Stokes  problem.  We  can  apply  the 
theory  of  the  mixed  finite  element  method  [2]  and  obtain  the  convergence  under  the  inf-sup  condition  ( 
Babuska-Brezzi  condition  )  similarly  to  the  case  in  the  Cartesian  coordinates.  Recently  Franca  et  al.  [3], [4] 
have  developed  the  stabilized  finite  element  method  where  the  inf-sup  condition  is  not  required  Any 
combination  of  velocity  approximating  space  and  pressure  approximating  space  is  possible.  We  show  that 
the  stabilized  method  can  also  be  applied  to  the  axisymmetric  problem  under  a  suitable  subdivision  of  the 
domain  and  a  proper  choice  of  base  fiinctions.  An  error  estimate  in  the  proper  norm  is  also  obtained. 
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NUMERICAL  SIMULATION  OF  FREE  BOUNDARY 
PROBLEMS 
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Introduction 

Boundary /Initial  value  problems  in  which  part  or  whole  of  the  flow  domain  is  not  known  apriori  and  needs 
to  be  determined  as  part  of  the  solution  procedure  are  called  free  boundary  problems.  Problems  of  this 
nature  are  of  interest  in  crystal  growth  processes,  open  channel  flows,  melting  and  solidification,  jets, 
flame  propagation,  groundwater  seepage,  metal  forming  processes,  draining  film  flows,  and  various  other 
fields  of  engineering.  Additional  non-linearity  is  introduced  into  an  already  non-linear  governing  equations 
due  to  the  unknown  free  boundary  position.  Due  to  this  reason,  analytical  treatment  of  these  problems  is 
extremely  difficult.  In  this  paper  a  numerical  model  developed  based  on  the  Arbitrary  Lagrangian 
Eulerian(ALE)  description  of  flow  is  described.  This  procedure  has  been  used  to  solve  open  channel  flow, 
thermocapillary  and  buoyancy  driven  flow,  and  draining  film  instabilities,  all  in  two  dimensions. 
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1.  INTRODUCTION 

An  upwind  finite  element  scheme  for  the  incompressible  viscous  flow  at  high  Reynolds  number  was 
proposed  by  Tabata  and  Pujima[l].  The  scheme  has  the  potential  to  appropriate  the  advection  term  in  third- 
order  accuracy.  We  apply  it  to  a  two-dimensional  non-stationaiy  analysis  of  air  flow  around  an  Automated 
Guided  Vehicle  (AGV).  The  results  are  at  least  qualitatively  good  and  compatible  with  experimental  ones. 
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1.  INTRODUCTION 

It  is  one  of  the  most  important  areas  in  technological  and  engineering  fields  to  calculate  transient 
flows,  such  as  those  found  in  Bows  of  soils,  cements,  metals,  and  polymers  etc.  There  are  so  called 
particle-in-cell  (PIC)  method(^)  markerand-ceil(MAC)  method(^)  analyzing  the  flows.  These  methods 
make  use  of  marker  particles  to  show  free  surfaces.  The  flow  patterns  can  be  visualized  if  the  markers  are 
also  arranged  inside  the  transient  flow.  These  methods,  however,  are  difficult  to  apply  under  complicated 
fixed  boundary  conditions.  Such  conditions  can  easily  be  covered  if  the  finite  element  method  is  used. 
Using  the  finite  element  method,  markers  can  be  also  used  for  the  analysis.  Shiojima  et  al.  have  introduced 
the  area  co-ordinate  system  of  the  linear  triangular  element,  which  is  utilized  for  finding  the  new  marker 
position,  into  a  four-node  rectangular  isoparametric  element.  The  two  kinds  of  interpolation  systems, 
however,  give  complicated  solution  algorithm  and  may  give  unexpected  errors  on  calculated  results.  Also, 
this  method  requires  liner  mesh  arrangement  to  obtain  an  accurate  result.  This  implies  that  a  large  amount 
of  extra  storage  may  be  required  for  the  analysis. 

in  this  study,  we  present  a  finite  element  method  using  six-node  triangular  isoparametric  element  for  the 
analysis,  in  which  marker  particles  are  introduced  to  represent  the  transient  creeping  fluid  flow  motion.  For 
determining  the  marker  positions  in  an  element,  area  co-ordinate  system  of  the  triangular  element  is  used. 
With  the  six-node  element,  the  determination  of  the  new  marker  position  in  an  element  becomes  very 
simple  compared  with  those  using  the  linear  triangular  element.  It  is  because  that  we  can  make  use  of  three 
mid-node  points  in  each  side  of  the  triangle  to  the  corresponding  three  area  co-ordinates.  Also  the  velocities 
at  the  six  nodes  obtained  by  the  calculation  can  be  directly  used  for  the  interpolation  of  the  marker 
velocity.  The  program  becomes  simple  and  higher  order  accuracy  can  be  expected,  and  course  mesh 
arrangement  can  be  available  with  the  present  method.  To  verify  the  scheme,  we  present  a  simple 
axisymmetric  problem.  The  time  dependent  creeping  flow  of  incompressible  Newtonjan  fluid  is  considered 
in  this  study.. 
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IDENTIFICATION  OF  THERMAL  CONDUCTIVITY  OF  GROUND  MATERIALS 
BY  NONLINEAR  LEAST  SQUARE  METHOD 
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I.  INTRODUCTION 

A  nonlinear  least  square  method  which  numerically  identifies  the  thermal  conductivity  of  ground 
materials  is  presented  in  this  paper.  The  performance  function  expressed  by  the  sum  of  squared  differences 
between  observed  and  computed  temperature  is  employed  to  estimate  the  thermal  conductivity.  The 
Fletcher  Reeves  method  is  chosen  as  the  methods  to  minimize  the  performance  function.  The  identification 
analysis  of  the  thermal  conductivity  of  the  ground  is  carried  out  by  using  the  temperature  measured  in  an 
actual  field. 
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1.  INTRODUCTION 

The  Crank-Nicolson(CN)  method  is  a  kind  of  useful  numerical  time  integration  methods  for 
simulating  some  transient  physical  phenomena  such  as  heat  or  mass  transfer,  wave  propagation  and  fluid 
flow  problems.  Although  the  method  has  some  advantages  as  unconditionally  stable  and  obtaining  second 
order  accuracy  in  time,  that  has  some  difficulties  to  cany  out  a  calculation  for  some  large  scale 
computations  since  that  is  a  kind  of  implicit  schemes.  The  author  and  co-workers  presented  the  block-wise 
local  CrankNicolson(LCN)  method  to  solve  one-  and  two-dimensional  heat  transfer  problems  with  the 
finite  difference  method[l].  The  method  inherits  some  advantages  from  the  CN  method  even  if  the  LCN 
method  is  formulated  with  an  explicit  form.  However  the  LCN  method  with  block-wise  decomposition  is 
not  extensible  to  three-dimensional  problems.  The  aim  of  the  paper  is  to  present  an  explicit  numerical 
integration,  which  is  unconditionally  stable  and  has  the  second  accuracy  in  time,  for  a  semi-discretized 
equation  which  is  derived  from  applying  the  finite  element  method  in  time  and  the  mass  lumping 
technique.  The  present  method  is  carried  out  with  explicit  calculations,  that  is,  the  method  is  formulated 
without  linear  solvers  as  Gauss-Jordan,  LU-decomposition,  Gauss-Sidel,  SOR  method,  etc..  Therefore  the 
author  call  the  present  method  the  explicit  Crank-Nicolson  method.  We  show  the  idea  of  the  present 
method  and  some  properties. 
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